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Abstract 

A calculational scheme for obtaining the electric polarizability of the neutron 
in lattice QCD with dynamical quarks is developed, using the background field 
approach. The scheme differs substantially from methods previously used in the 
quenched approximation, the physical reason being that the QCD ensemble is no 
longer independent of the external electromagnetic field in the dynamical quark 
case. One is led to compute (certain integrals over) four-point functions. Par- 
ticular emphasis is also placed on the physical role of constant external gauge 
fields on a finite lattice; the presence of these fields complicates the extraction of 
polarizabilities, since it gives rise to an additional shift of the neutron mass unre- 
lated to polarizability effects. The method is tested on a SU (3) flavor-symmetric 
ensemble furnished by the MILC Collaboration, corresponding to a pion mass of 
to^ = 759 MeV. Disconnected diagrams are evaluated using stochastic estimation. 
A small negative electric polarizability of a = (—2.0 ± 0.9) • 10 -4 fm 3 is found for 
the neutron at this rather large pion mass; this result does not seem implausible 
in view of the qualitative behavior of a as a function of m n suggested by Chiral 
Effective Theory. 
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1 Introduction 



An important characteristic of hadrons is their stiffness when subjected to outside forces, 
which are typically conveyed by external electromagnetic fields. This response is sum- 
marized in hadron polarizabilities. Understanding these quantities will contribute to 
making hadron structure more palpable. Experimentally, polarizabilities are accessible, 
e.g., via soft Compton scattering; heuristically, in such an experiment, the photon elec- 
tric and magnetic fields polarize the target hadron, which in turn manifests itself in 
the Compton scattering amplitude observed. Accordingly, polarizabilities are effects of 
second order in the external fields. 

The aforementioned sensitivity of low-energy Compton scattering to hadron structure 
can be cast in precise terminology [1], permitting stringent tests of theoretical under- 
standing of that structure. Starting with the leading order in the low-energy expansion, 
the non-Born (i.e., structure-dependent) part of the scattering amplitude is determined 
by the static dipole electric and magnetic polarizabilities a and (3. These are given by 
the hadron mass shift in the presence of external static electric and magnetic fields, 
specifically the part of the mass shift which depends quadratically on those fields, in 
accordance with a (spin-independent) effective dipole interaction Hamiltonian 



The present investigation focuses on the electric polarizability a of the neutron. 

Lattice hadron polarizability calculations have hitherto been carried out only in the 
quenched approximation [2-9]. The reason for this lies in the fact that, in the case of 
polarizabilities, the complication implied by going from a quenched to an unquenched 
calculation involves more than just the usual vastly increased effort required to generate 
a dynamical quark ensemble. In addition, a quenched calculation is simpler due to the 
gauge ensemble being independent of the external electromagnetic field; after all, the 
only way the external field can influence the gauge ensemble is through the quarks, 
whose backreaction on the gauge fields is precisely truncated in a quenched calculation. 
This is no longer true in the dynamical quark case. 

This physical difference manifests itself formally in the fact that substantially differ- 
ent computational schemes have to be used in the dynamical quark case as compared to 
the quenched case. In the quenched case, one can simply generate gauge configurations 
in the absence of the external electromagnetic field and introduce the latter a posteriori 
by an appropriate modification of the link variables in those configurations. The req- 
uisite hadron two-point functions are then evaluated directly using the modified gauge 
configurations. By contrast, in a fully dynamical calculation, as discussed in more detail 
below, one in principle would need to generate the gauge ensemble anew for each ex- 
ternal field considered. The prohibitive cost of such a scheme can be mitigated to some 
extent by expanding in the external field, leading, in effect, to the calculation of (certain 
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space-time integrals over) four-point functions. In general, these include disconnected 
contributions. However, even resorting to such a four-point function method leads to 
a substantially more expensive calculation than one is confronted with in the quenched 
case; an early exploratory study of four-point function methods [3] (using a quenched 
ensemble) highlights this poinlo 

The most of this situation has been made hitherto in a series of investigations [2,5-9] 
taking full advantage of the simplifications offered by the quenched approximation. An 
initial study of the electric polarizability of neutral hadrons [2] using staggered fermions 
yielded results both for the neutral pion and the neutron. This was later extended to 
include all neutral members of the baryon octet, as well as the baryon decuplet and the 
vector meson octet, in investigations using both Wilson and clover fermions [5,7]; these 
studies surveyed a range of pion masses down to about 500 MeV. Also the magnetic 
polarizability of a wide array of hadrons was investigated using the same range of pion 
masses and fermion actions [6,8,9]; these studies included also charged hadrons, in 
particular the entire baryon octet and decuplet as well as selected pseudoscalar and 
vector mesons. 

The present work extends the aforementioned work in two main aspects: 

• Use of a dynamical quark ensemble: As discussed above, dynamical quark calcu- 
lations of polarizabilities were rendered intractable in the past by the associated 
computational cost. Recent increases in available computing resources are making 
quantitative four-point function calculations, appropriate for dynamical ensembles, 
feasible. This investigation presents the first result for the electric polarizability of 
the neutron in a dynamical quark ensemble, albeit obtained at a still rather heavy 
pion mass of 759 MeV. 

• Recognition that, on a finite spatial volume, a constant gauge field is not a pure 
gauge, but has physical consequences which must be disentangled from polarizabil- 



Some elaboration on the latter issue, which is also relevant in the quenched approxi- 
mation, is useful at this point. As noted further above, hadron polarizabilities can be 
probed via the mass shift in the presence of external electromagnetic fields. The hitherto 
preferred method [2, 5, 7] of introducing a constant electric field in, say, the 3-direction 
is to represent it by a non-vanishing 3-component of the gauge field, 



This choice has the advantage that jumps in the gauge field at the lattice boundaries 
(inducing spurious localized electric fields there) occur only in the temporal direction. In 

lr The four-point functions considered in [3] are different from the ones which are calculated in the 
present work, since different theoretical approaches are used. However, the computational complexity 
resulting from the two approaches is similar. 



ity effects. 
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this case, hadron two-point functions evaluated in the bulk of the lattice are insensitive 
to the jumps, which only occur far in the past or the future from the point of view of 
the measurement. 

However, there is an ambiguity in the prescription ([2]), namely, at which time to 
one chooses to begin counting time. Different choices of t correspond to different 
constant shifts of A 3 . Working in a spatially infinite setting, this ambiguity would 
be inconsequential, since constant gauge fields are then pure gauges. However, on a 
finite space, the spatial boundary conditions (which in the following will be taken to be 
periodic) restrict the available gauge transformations and only allow for discrete shifts of 
the gauge fields. As a simple example, consider a charged particle in a constant field on a 
circle of length L described by the Hamiltonian H = (— id x + A) 2 . Its energy eigenvalues 
are E n = (2irn/L + A) 2 , where n can be any integer. The ground state energy therefore 
is E = A 2 as long as A e [— ir/L, tt/L], and E is periodic in A with period 2n/L, 
reflecting the residual discrete gauge invariance. The spectrum explicitly depends on A 
and, compared with the case A = 0, the ground state energy can deviate by as much as 
AE = tt 2 /L 2 . Thus, while this is ultimately nothing but a finite size effect, it vanishes 
rather slowly (only as a power of L) as the spatial volume is increased. It is a priori 
unclear how difficult it is in practice to deal with this effect by using different lattice 
sizes. Within the present investigation, that avenue is closed at any rate, since the 
dynamical quark ensemble which will be used is only available at one spatial volume. 
Instead, measurements at several different t in ([2]) will be used in order to treat this 
effect. 

Another aspect of the same issue is that the Hamiltonian in the presence of the field 
(T5]) is not time-independent. Physics at two widely separated times differ precisely by 
a shift in the external gauge field A3. If the electric field E is very small, the strong 
dynamics can instantaneously adjust to the change in the external field as time passes; 
one will observe an adiabatic change in the physical spectrum. The hadronic two-point 
function will generally not fall off as a simple exponential in Euclidean time, even for 
large such times. Both the hadron's energy as well as its wave function will contain time 
dependences. This is reminiscent of the behavior induced by the acceleration of charged 
hadrons in the electric field [10]. That particular effect is negligible for sufficiently heavy 
hadrons, such as nucleons. By contrast, the constant gauge field effect discussed here 
is one of the dominant effects, and care needs to be taken to disentangle it from the 
nucleon polarizability. 

A way to avoid the time dependence discussed above is to instead use a gauge field 
representation of the type 

A = -Ex 3 (3) 

inducing the same external electric field as ([2]); indeed, in the present work, also this case 
will be investigated. However, it should be noted that this choice also has disadvantages 
which, on balance, are no less problematic than the time dependence engendered by ([2]). 
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Namely, the advantage of time-independence using ([3]) is offset by the fact that spatial 
translational invariance is lost; the representation ([3]) conflicts with spatial periodicity 
and the periodic boundary conditions enforce a spike in the electric field resulting from 
the jump in A as one crosses the boundary of the lattice in the 3-directiorj§. In effect, 
the neutron thus propagates in a spatially varying potential and its energy contains, 
e.g., contributions from quantum mechanical zero-point motion in that potential. Also 
the neutron's internal wave function is distorted by the electric field spikes. As a con- 
sequence, it is not straightforward to isolate the polarizability from the full measured 
mass shift. 



2 Measurement method 
2.1 Neutron two-point function 

The objective of the present investigation is to extract the neutron mass from the neutron 
two-point function in the presence of an external electric field. The neutron two-point 
function is the correlator 

(N a ,(x')N a (x)) = \ [ [DU}[D^}[D^}exp(-S[^^,U})N a/ (x')N a (x) , (4) 



Z 

with the lattice discretization of the functional integral to be specified below. Both 
the action S and the (smeared) neutron fields N, N in general depend on the external 
electromagnetic field A^. At face value, this would imply that one needs to generate 
lattice ensembles using an action modified by the external field in order to evaluate 
03]). That would clearly be prohibitively expensive using dynamical quarks. However, 
decomposing the action as 

S = S + S E , (5) 
where So is the action in the case of vanishing external field, one can rewrite (JH) as 

<*,<*-)*.<«)> = <£f!W^ (6) 

(e bE )o 

with (. . .)o denoting the average in the absence of the external field, 

(O) = 1-J [DU][D$\[Di/>] exp(-So)0 . (7) 



2 Note that, in the lattice formulation, there exist discrete choices of E which mitigate this problem, 
namely, integer multiples of 2n/aL, where L is the extent of the lattice in the relevant direction and 
a denotes the lattice spacing. However, this argument relies on the compactness of the gauge link 
variables and, on realistic lattices, corresponds to strong electric fields. By contrast, hadron electric 
polarizabilities are given specifically by the term quadratic in E of a Taylor expansion of their mass. To 
isolate this term, it is necessary to vary E over a denser set of values than provided by the aforementioned 
discrete choices, for which the Taylor expansion will generally not converge well on lattices of a practical 
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While this reduces the problem to integrations over the lattice ensemble in the absence 
of the external field, expectation values such as the ones in generally suffer from 
severe overlap problems. The crucial step which renders the problem somewhat more 
tractable results from the fact that it is sufficient to know the quadratic term in the 
Taylor expansion of ([6]) with respect to the external field in order to extract the neutron 
polarizability. Then, one can expand 



and the evaluation of ([HD reduces to the calculation of certain space-time integrals over 
four-point functions, as will be discussed in detail below. Before proceeding to describe 
this perturbative expansion, it is now possible to specify how the functional integration 
in ((7j) will be carried out. As usual, decomposing Sq into its pure gauge and fermion 
parts, 



the integration is cast in terms of an average over an ensemble of gauge fields U, 



governed by the action Sa,eff which includes the effects of both the pure gauge term as 
well as the determinant of the Dirac operator from Sp. Adopting 2 + 1 flavor Asqtad 
quark fields to evaluate the determinant, one can utilize the corresponding dynamical 
quark ensembles made available by the MILC Collaboration [11-13]. The numerical re- 
sults reported in the present work were obtained using 99 configurations from the SU(3) 
flavor-symmetric ensemble with quark masses given by am s = ami = 0.05, where the lat- 
tice spacing a = 0.124fm is determined by heavy quark spectroscopy [14]. Computation 
at such a relatively large quark mass is comparatively inexpensive and serves primarily 
to validate the concepts developed in this work as well as giving a first indication of the 
feasibility of a broader calculational effort within the framework advanced here. 

The aforementioned configurations were originally generated on 20 3 x 64 lattices. 
In the present work, these lattices were chopped in half in the time direction, i.e., 
measurements were carried out on 20 3 x 32 lattices containing the first 32 time slices 
of the original 20 3 x 64 lattices. Correspondingly, quark propagators determining the 
quantities (O)u, cf. (TTTT) . were evaluated using Dirichlet boundary conditions at the 
temporal edges of the chopped lattices. The lattices were furthermore HYP-smeared [15] 
to reduce the effect of dislocations. 

For a given gauge configuration U, the expectation value 



exp(-S E ) = 1-S E + S 2 E /2 + ... 




So — Sg + Sf , 



(9) 




(10) 



(O)u 



J[DnDjj]exp(-S F )0 
f[D$][Dij>] ex P (-S F ) 



(11) 
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needs to be evaluated. At this point, a hybrid approach is adopted [16, 17]: While 
the gauge ensemble used in the average ffTUj) is generated using Asqtad quarks, (fTTj) is 
evaluated using domain- wall quarks [18,19]. The reason for this choice lies in the longer- 
term goal of extending the present investigation to light quark masses at which a chirally 
well-behaved quark discretization becomes important. Thus, the fields ip, ip in ffTTl) are 
taken to be five-dimensional, ip(x, s) and ip(x, s), where x labels four-dimensional space- 
time and s labels the fifth coordinate. The latter is subdivided into L s = 16 spacings, 
s G {0, . . . , L s — 1}, and ip(x, s) = outside that range (i.e., derivatives in the fifth 
direction have hard boundaries). This value of L s is sufficient to keep the residual mass 
m res which characterizes explicit chiral symmetry breaking by the domain-wall fermion 
discretization suppressed by more than an order of magnitude compared to the quark 
mass discussed further below [16, 17]. The boundaries s = and s = L s — 1 provide 
domain walls which support quasi-f our- dimensional light physical quark modes; left- 
handed modes are bound to s = and right-handed modes to s — L s — 1. It is useful 
to define corresponding four-dimensional projected quark fields 



= ^- Y ^ij(x,0) + ^^ij(x,L s -l) 
*(x) = i>(x,0)^^ + ^(x,L s -l)^^ 
In terms of the above fields, the action Sp in ffTTl) reads 



(12) 
(13) 



X,S fl 



\U^{x)il){x + e M , s) - ip(x, s)) 



-LfUx - e^)i)(x - e M , s) + ip(x, s) 



^{x, s)M 5 ip(x, s) + ^2 ty(x)mfty( 



x 



(14) 



where \x runs over all five dimensions and U^ = l. Note that the fermion fields also carry 
a flavor index; in the final term, i.e., the quark mass term, which is constructed directly 
in terms of the projected four-dimensional quark fields ^ and rrif in general represents 
a (diagonal) matrix in flavor space. In the SU(3) flavor-symmetric case studied here, 
rrif is given by one single number. 

Using spectral flow analyses, the five- dimensional mass parameter M 5 in ffTTl) was 
chosen to take the value M5 = 1.7 in order to optimize the chiral properties [16,17]. 
Finally, the quark mass was adjusted such as to match the pion mass obtained in the 
present hybrid approach to the lightest pion mass extracted from a pure Asqtad calcu- 
lation [14]; this yields [16,17] the choice arrif = 0.081. 

The domain wall fermion action ffTTl) also determines the interaction between the 
quarks and the external electric field. The additional electromagnetic gauge field 
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generating the external electric field modifies the gauge link variables, 

U^x) — > ex^{iaq f A^(x))U^{x) , (15) 

where a denotes the lattice spacing; note that the fractional electric charge qj varies ac- 
cording to flavor. Note also that the particular forms of used in this work, cf. (J2J), ([3]), 
are all such that A^ is constant in the /i-direction; hence the simple form f|T5|) for the 
exponentiated integral along the link. Inserting the modified link variables (j!5p into the 
domain wall fermion action ( fT4l) and separating off the part which remains for vanishing 
external field, A^ = 0, yields the five- dimensional interaction 

Se M = -EE^^)(^ i (e w ^ (x) -l)^(^(x + e^, S ) (16) 

X,S fl \ £ 

+ l A^ L (e~^f A ^-^ - l) Ul(x - e^{x - e„ s)) 

generating a vertex which couples the five-dimensional domain wall fermion fields if), ip 
to the external field. 

To arrive at a practicable computational scheme, in the calculations presented further 
below, the external gauge field A^ is not coupled directly to the five-dimensional fields 
according to ffTBl) . but instead to the corresponding four-dimensional projected quark 
fields Accordingly, a renormalization factor zy must be included with the four- 

dimensional coupling to compensate for the effect of the projection of the quark fields. 
Thus, the modified interaction vertex used in practice is 

Se = -zvEE*^)(^^( eia9MM(a:) - 1 )^(^)*(^ + e M) (17) 

+ i±3i (e-^fM*-*.) _ i) Ufa _ e ^( x _ e ^ . 

The renormalization factor zy will be determined in section 14.11 The reason for the 
adoption of the modified interaction ([17]) lies in the practical expense of storing full five- 
dimensional propagators as opposed to ones which have been projected down to four 
dimensions at source and sink. This modus operandi constitutes a compromise which 
certainly should be revisited as storage constraints change. Using full five-dimensional 
propagators and coupling the conserved five-dimensional current to the external elec- 
tromagnetic field directly via f|T6l) would be the most consistent treatment, and would 
eliminate the need for renormalization of the interaction vertex. 

Finally, it is necessary to specify the neutron sources and sinks N, N in (j3j): 

N a (x) = 6 aP (C l5 ), s e bcd Q^(x)Q^(x)Q^(x) (18) 
N a (x) = Q { £{x)Q%\x)Q${x)e M 6 a p(C 75 ) lS , (19) 
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where C denotes the charge conjugation operator and Q is a Wuppertal-smeared [20] 
quark field (with the superscript denoting flavor), constructed iteratively as (where the 
superscript now momentarily labels iterations): 

Q«(x) = (1 - G^Q^ix) + a £ V\{x - e^Q^x - e M ) . (20) 

Here, a and the number of iterations i max are free parameters, chosen such as to generate 
a good overlap between the neutron source and the true neutron ground state [21]. 
The iteration starts at = ^ and ends at Q^ lmax ) = Q. The sum over directions 
/i in (1211 runs only over the three spatial dimensions, but includes terms associated 
with both positive and negative displacements in each dimension (i.e., e-^ = — e M , 
Ul. ix {x — e_ M ) = U^(x)). Note that smearing constitutes a linear operation on the quark 
fields, i.e., there exists a matrix P such that 

Q(x) = P(x, y)*(y) , Q(x) = *{y)P\y, x) . (21) 

P is proportional to the unit matrix in the Dirac indices, but not in the space-time and 
color indices, nor in the flavor indices once the external electric field is introduced via the 
substitution ffT5]) . Note, thus, that the presence of the external electric field can influence 
the smearing if one insists on manifest invariance of the neutron sources and sinks with 
respect to gauge transformations of the external field. However, it is not imperative to 
preserve such manifest invariance; not doing so merely corresponds to evaluating (gauge- 
invariant) physical observables in a particular gauge. In the treatment to follow, the 
most general case will be considered, i.e., the perturbative expansion discussed below 
will include the diagrams resulting from expanding the source and sink fields in the 
external field. This will make it possible to separately assess the influence of such 
terms. Ultimately, unambiguous extraction of the neutron electric polarizability will be 
seen to necessitate discarding such diagrams, and thus foregoing manifest invariance of 
the neutron sources and sinks with respect to gauge transformations of the external field; 
nevertheless, it will be verified that the effect of including additional smearing diagram 
contributions on the final result for the polarizability is negligible, thus rendering this 
issue moot in any case. 

2.2 Perturbative expansion 

Having defined all of the objects entering the neutron two-point function (j3J), one can 
proceed to extract the quadratic term of its Taylor expansion with respect to the external 
field A p. Both the interaction Se and the smeared neutron sources N and N in general 
contain a dependence on A^. Expanding (TTT|) . one obtains two relevant vertices, 

S e = Se,i + S e ,2 + o(A 3 ^ , (22) 
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which can be written as bilinear forms, 

3 E)i = *M<* , (23) 

with 

Mi(x,y) = -iaz v q f J2 f 1 ltl A lJL (x)U lJL (x)5(x + e^y) (24) 



1 + 7/, 



2 

,2 



4* fa - ejUfa - e fl )5(x - e M , y)j 



M 2 (x,y) = ^zyq^i^Y^A^U^Six + e^y) (25) 



^(^-e M )^(a;-e M )5(x-e M ,y)^ . 



Thus, Mi and M 2 are matrices in the space-time, color, Dirac and flavor indices, sum- 
mation over which is implied in fl2"31 . 

On the other hand, also the smeared fields defined by (120]) need to be expanded in 
the external field, 

Q (l) = Qf + Qf + Qf + o Ul) (26) 



(where the subscript denotes the order in the external field). Modifying the link variables 
in (|20|) according to (|T5|) and expanding in A^, one has an iterative construction of the 
smeared fields separated order by order in the external gauge field, 

Qf{x) = (l-6a)Qt 1 \x)+ajrul(x-e fl )Qt 1) (x-e fl ) (27) 

Q?(x) = (1 - 6a)Qt 1] (x) +ajr U^x - e M ) fat^x - e») (28) 

fi=±i 



-iaqfA^x - e^Qfr l \x - e A 

Qf{x) = (l-6a)Qt l \x)+af:ul(x-e,)(Qt 1 \x-e») (29) 

ii=±i 



-iaqjA^x - e^)QX l \x - e 



a 2 q) 



2 A^ix — e M )Qo \x — e M )J . 

Equivalently, the smearing matrix P in (I2T]) can be written in expanded fashion, 

P = P + P x + P 2 + O (Al) (30) 
(the original quark field is of course of zeroth order in A^). 
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Returning to the neutron two-point function, expanding (jBD in powers of Se, inserting 
(j22p and discarding terms which contribute only at higher than quadratic order in the 
external field yields 

(N a ,(x')N a (x)) = ((l-SE + Sl/2)N a/ (x')N a (x)) o (l + (SE-S 2 E /2) + (S E )l) 
= (N a ,(x')N a {x)) 

[S E ,i + S E ,2 - S 2 E J2) N a ,(x')N a (x)) o 
[Sb,i + S E ,2 - S 2 EA /2) o (N a ,(x')N a (x)) n (31) 
-(S E ,i) (S E ,iN a ,(x')N a (x)] 
+ {S E A)l(N a >(x')N a (x)) o . 

As usual, the denominator in the original expression (jH]) has the effect of subtracting 
disconnected (in the statistical sense) pieces. Furthermore, inserting the more specific 
forms (fI8]l , (|T9|) , (l2~Tl) and (|23|) , one arrives at (the superscripts of the smearing matrices 
P and the quark fields \1/ denoting a fixed flavor): 

(N a >(x')N a (x)) = S al p>(C^ 5 )^s'e b ' c ' d 'e bcd S al 3(Cj B ) 7S 

x*. u>)P¥?,{x>, v>)pM,(x>, w')P^\w, *)P/ C ( %, x)PH d \u, x) 
x - (^a(^0^^^^0^L%(^0^1(^)^i 7 d V)^g(«) (32) 
x f(f Mi*) + (f M 2 *) - -(f M!*)(f Mi*)^\ 

- (*%(u')^,(v>)¥^ (M^)) n 

x ((M,$) 



x + ^(*Mi*) + (*M 2 *) - -(*Mi*)(*Mi*)^ 

+ (($ M x *))* 

Applying Wick's theorem (i.e., evaluating the (. . .)u averages over the quark fields, 
cf. (fT0|) .(flTl)). and retaining only contributions quadratic in the external gauge field, 
one finally arrives at a diagrammatic representation, depicted in Fig. [TJ for the desired 
quantity, namely, the quadratic term in the Taylor expansion of the neutron two-point 
function with respect to the external field. The diagrams in Fig. [1] are to be read as 
follows: 

a. Solid lines are point-to-point quark propagators K^ c c (x f , x) = fif^h (x')^f^ {x))u- 
Note that these are propagators between four-dimensional sources and sinks, i.e., an 
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O O 




Figure 1: Contributions to the neutron two-point function quadratic in the external 
gauge field. The nomenclature is explained in detail in the main text. 
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initially four-dimensional source is propagated in five dimensions using the domain 
wall quark action (TH1) and finally projected back to four-dimensional space-time 
according to the correspondence f|T3|) . Of the three propagator chains connecting 
neutron source and sink, two are associated with flavor down, / = d, and one with 
flavor up, / = u. The quark loops imply a sum over all three flavors. 

b. O denotes quark source or sink smearing of zeroth order in the external field, i.e. 
Pq in the decomposition (I5U1) . Similarly, corresponds to Pi and # corresponds 
to P 2 . 

c. X denotes a vertex insertion linear in the external field, i.e., multiplication by M\, 
cf. ^2M . Similarly, ® corresponds to M 2 , cf. (125]) . 

d. At the neutron source and sink, symbolized by the ovals, color and Dirac indices 
must be contracted in accordance with the first line of the right-hand side of (I32p . 

e. Each diagram summarizes several elementary terms in the Wick expansion of fl32l . 
For every contribution in which the quark lines run literally as shown, there is a 
corresponding contribution in which the two down quark lines connecting neutron 
source and sink cross (i.e., the sinks are exchanged). The latter contribution re- 
ceives an additional minus sign from the exchange. Furthermore, for each diagram, 
there are several ways of distributing vertices and smearings over the quark sources, 
sinks and propagators, only one of which is shown in each case. Note that the com- 
binatorics are different for vertices and smearings. On the one hand, there are six 
ways of distributing two X vertices such that they reside on different quark lines 
connecting neutron source and sink (i.e., exchanging two such vertices amounts to 
a new contribution!); on the other hand, there are only three ways of distribut- 
ing two sink smearings in the neutron sink (since the product (Pq + Pi + P2) 3 
contains only three terms consisting of two factors Pi and one factor Po), and 
analogously for the neutron source. Note that the labeling of the diagrams reflects 
these multiplicities; in each label, the integer inside the parentheses denotes the 
number of individual contributions from the Wick expansion of (E?2"|) summarized 
by the diagram. 

f. As usual, each quark loop implies an additional minus sign. In order to keep with 
standard nomenclature, these signs were not absorbed into the prefactors, but must 
be included separately when evaluating the diagrams. 

g. Each contribution finally must be averaged over the gauge ensemble, where, as 
already remarked after eq. (T5TT) . statistically disconnected parts are subtracted. 

3 To be completely precise, this only applies when both vertices reside on quark lines connecting 
neutron source and sink; on the other hand, in the diagrams labeled J01(2), J02(6) and -£T01(2), no 
additional contributions stemming from exchange of the vertices are implied. Any such duplications 
which may arise from the Wick expansion of (|32|) are already taken into account through the prefactor 
of the diagram (such a duplication actually only occurs in the case of J02(6), the statistically connected 
part of which, cf. item g., originally enters with a prefactor of 1/2). 
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Thus, denoting the gauge ensemble average as 

(0)g = ^r f \DU] exp(-£ G>e// )0 , (33) 
diagram J01(2) is to be evaluated as 




(where of course only items a.-f. apply to the objects inside the averages). Diagrams 
J02(6), ,703(2), Jll(6) and J12(6) are treated analogously The more complicated 
case K01(2) is evaluated as 
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Note that the naming of the different contributions is intended to be mnemonic. The 
initial character differentiates between contributions of varying number of disconnected 
parts; I denotes connected diagrams, J disconnected ones with two parts and K the 
disconnected diagram with three parts. The digit following the initial character indi- 
cates the power of the external electric field contributed specifically by the source and 
sink smearings. The next digit is simply a running index numbering the contributions 
in each class. Finally, as already mentioned under item e. above, the integer in the 
parentheses denotes the number of individual contributions from the Wick expansion of 
(|32l) summarized by the diagram. 
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2.3 Calculational details 



The code written to compute the diverse diagrams in Fig. [T]relied heavily on the Chroma 
Library for Lattice Field Theory [22,23]. In practice, the propagator chains connecting 
neutron source and sink were calculated in sequential fashion. Starting from a space- 
time location x and a specific set of color, Dirac and flavor indices a, a and /, as well 
as choosing the desired order % of the smearing in the external field, one constructs 
the smeared source vector 5p a P- b a (y, x) using the appropriate iterative procedure (I27|) . 
( |28|) or (|29|) . While no loss of generality is incurred by performing the calculation 
for only one particular x, all combinations of the other indices are ultimately required 
for the contractions at the neutron source (of course, different flavors are related in a 
trivial manner). Propagating the aforementioned specific smeared source vector yields 
directly the smeared-to-point propagator K^J cb (z,y)P-^ a (y,x). A vertex insertion im- 
plies multiplication with the corresponding matrix Mj, yielding a new source vector 
Mj^^(w, z)K^Q cb (z,y)P-^ a (y, x). This source vector is then again propagated^, thus 
building up the propagator chain sequentially. When finally arriving at the neutron 
sink, the appropriate sink smearing is applied, using again ( I2T1) . (I2"gj) or (|2"§1) . 

The disconnected quark loops were evaluated using stochastic estimation. To esti- 
mate the trace over all indices implied by the loop, a basis of 120 stochastic sources 
(240 for two cases of external fields which engender particularly strong statistical fluc- 
tuations, cf. section I4T21 was used. Again, starting at each stochastic source, propagator 
chains were constructed sequentially, and finally contracted again with the stochastic 
source. The sources were complex Z(2) sources, distributed homogeneously over space- 
time, Dirac and color spaced, i.e., each point in that product space was associated with 
a value from the set 

{l+z,l-z,-l + z,-l-z} (34) 

with equal probability. 

Disconnected contributions exhibit strong statistical fluctuations, and two possibili- 
ties of reducing these fluctuations were investigated. On the one hand, the consequences 
of only switching on the external electric field a short time before the introduction of 
the neutron source and switching it off soon after the annihilation by the neutron sink 

4 Note that the positions of the interaction vertices in the diagrams in Fig. [T] are not external parame- 
ters, but integration variables. Thus, e.g., diagram 701(12) does not represent a full four-point function, 
but only a very specific space-time integral over a four-point function. It is these integrations which 
render the calculation tractable by the sequential procedure described here; they provide precisely the 
contraction between a vertex and an attached propagator which permits treating an inserted vertex 
simply as one single new source, devoid of external parameters and spread out over all of space-time, 
to be submitted to the subsequent propagation. 

5 For the flavor SU (3)-symmetric ensemble m u — m c i = m s used in this work, it is sufficient to 
consider one flavor and weight the result by the appropriate combination of fractional charges to obtain 
the full value of the loop diagram. 
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Figure 2: Comparison of results obtained using chopped and unchopped external fields, 
as described in the main text. Left panel displays the sum of the connected diagrams 101, 
102 and J03; right panel the sum of the disconnected diagrams J01 and J03. Results 
are shown as a function of temporal source-sink separation, in each case normalized by 
the neutron two-point function in the absence of the external field, i.e., shown are the 
contributions by the respective subsets of diagrams to the ratio defined in eq. ( 140]) . 
All measurements are taken at integer times; data are slightly displaced from those 
times in the figures for better readability. The electric field E providing the scale is 
cast in Gaufiian units. Shown are unrenormalized raw data, i.e., for the purpose of this 
comparison, z v = 1 in the vertices flUJ) . fl25|) . 



were explored. This procedure will be referred to as "chopping" the external field in 
the following. It is motivated by the expectation that, if sufficient time has elapsed 
between the introduction of the neutron source and the neutron mass measurement to 
filter out the true neutron ground state, then also any switching-on effects generated 
prior to the introduction of the neutron source will have decayed. However, the sta- 
tistical fluctuations of disconnected diagrams will be significantly affected by chopping 
the external field. Summing up contributions due to the coupling of the external field 
to vacuum fluctuations far in the past or the future of the neutron mass measurement, 
while not expected to influence the outcome of the latter, will certainly add statistical 
noise to it. Chopping the external field can reduce that noise significantly by discarding 
irrelevant vacuum fluctuations. This was tested using the external field A 3 = E(t — t ) 
with to = —10a, where here and in the following, the temporal lattice boundaries are 
located at t = —10a and t = 22a, and the neutron source is located at t = 0. Figs. [2] 
and [3] compare results obtained without chopping, i.e., A 3 = E(t — t ) throughout the 
lattice, with results obtained by setting A3 = E(t — to) only for —a < t < 14a and 
A3 = for other times. As expected, no significant differences arise in the measured 
correlator ratio. This is particularly clear in the connected contributions, which are 
determined very accurately; in the disconnected contributions, a significant reduction 
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Figure 3: Comparison of results obtained using chopped and unchopped external fields, 
analogous to Fig. but showing the contributions of the connected diagrams 701, 102, 
103 and the disconnected diagrams J01, J03 in one figure. 



of the statistical uncertainty results. Note that the neutron mass shift is ultimately 
extracted specifically from the slope of the correlator ratio shown in Figs. [2] and [3j as 
discussed in section |3j Due to the advantages offered by chopping the external field, 
all further measurements reported in the following, cf. in particular section 14.21 were 
obtained using chopped external electric fields. 

On the other hand, a further possibility of reducing the uncertainty of stochastic 
estimation which was explored is dilution [24], specifically dilution in the Dirac index. 
In other words, besides the stochastic estimation scheme described above, also an al- 
ternative scheme was considered in which each value of the Dirac index in the loop 
trace was considered separately, with Z(2) sources distributed homogeneously only over 
space-time and color space in each case, the sum over those values yielding the Dirac 
trace at the end. The comparison between the two schemes was carried out for the 
external field A 3 = Et, chopped as described further above, with the neutron source 
again located at t = 0. Fig. 0] shows the respective results obtained for the disconnected 
diagrams J01 and J03. Evidently, for this particular external field and dilution scheme, 
there is no computational advantage in dilution; the statistical uncertainty in fact is 
slightly larger in the diluted case. As a consequence, dilution was not considered any 
further in the present investigation; this does not exclude that a comprehensive survey 
of various implementations of dilution could yield computationally more advantageous 
schemes. 



3 Interpretation of the neutron two-point function 

The standard method of extracting ground-state hadron masses is to project the hadron 
two-point function onto a definite momentum, consider an appropriate Dirac component, 
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Figure 4: Comparison of results obtained using undiluted and diluted stochastic sources, 
as described in the main text. Diagram J01 (left) and diagram J03 (right) are displayed 
as a function of temporal source-sink separation, each normalized by the neutron two- 
point function in the absence of the external field, i.e., shown are the individual con- 
tributions by the two diagrams to the ratio R 2 defined in eq. ( )40|) . All measurements 
are taken at integer times; data are slightly displaced from those times in the figures 
for better readability. The electric field E providing the scale is cast in GauBian units. 
Shown are unrenormalized raw data, i.e., for the purpose of this comparison, zy = 1 in 
the vertices (El), (BSD- 



and compare the measured data to the corresponding spectral representation. Choosing, 
specifically, zero momentum and unpolarized neutron states, 

G{p = 0,t)= J dVTr (^^-(N(x')N(x))) — ► W exp(-mt) (35) 

for sufficiently large times t, where the neutron source location defines t = 0, and W 
characterizes the overlap between the state created by the operator N and the true 
neutron ground state. Thus, the neutron mass m can be extracted from the exponential 
decay of the correlator (1351) . 

Furthermore, if one is calculating the correlator as a function of a small external 
parameter, such as an external electric field E, one can expand in E, 

m = m + m x E + m 2 E 2 + O (e 3 ) (36) 
W = W + W X E + W 2 E 2 + O (E 3 ) (37) 



and then the Taylor expansion of (135]) contains the quadratic term 

G( 2 )(p = 0,t)= exp(-m t) [W 2 - W m 2 t - W x m x t + W m\t 2 /2\ E 2 . (3* 
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Using the fact that the neutron's electric dipole moment vanishes, mi = 0, and dividing 
by the correlator Go obtained in the absence of the external field, 

G (p = 0, t) — > Wo exp(-m t) , (39) 

one has 

G^(p = 0,t) fW 2 \ 2 

allowing one to extract the neutron electric polarizability 

a = -2m 2 , (41) 

cf. (0Q), from the slope of (|4"0|) as a function of t. 

Two assumptions underlie this procedure, namely, time-independence of the Hamil- 
tonian and spatial translational invariance. As already indicated in section (TJ neither of 
the external gauge fields (j2J) and ([3]) investigated in the present work satisfies both of 
these assumptions simultaneously. As a result, the standard analysis discussed above, 
which would be appropriate in infinitely extended space-time, needs to be reconsidered 
in more detail. 



3.1 Temporally varying gauge field 

Consider first the case of the external field ([21) , 

A 3 = E(t - to)=A + Et. (42) 

In this case, one does have spatial translational invariance, but there is no invariance 
under arbitrary temporal shifts. A translation in time corresponds to a shift in the 
constant component A of the gauge field, and, on a space of finite extent, different A 
are in general physically inequivalent, since only gauge transformations which shift A 
by certain finite increments exist. Therefore, the Hamiltonian in the presence of the 
field (1421) is time-dependent (with a periodicity which depends on the strength of the 
electric field E). 

These observations affect the above analysis in two ways. For one, the correlator (1351) 
really depends on two external parameters, E and A. As a consequence, the correlator 
ratio (j4"U|) contains all quadratic dependences which can result in the presence of the two 
parameters, 

R 2 (t) — _L (w 2 AA A 2 + W 2 AE AE + W 2 EE E 2 ) - (m AA A 2 + m AE AE + m EE E 2 ) t . 

(43) 

6 Note that one cannot argue external fields of the type ((2|) to be gauge-equivalent to time- independent 
ones such as ([3]) on a finite coordinate space; the corresponding gauge transformation conflicts with the 
boundary conditions. 
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On the other hand, for a small external field, the ground-state neutron mass and wave 
function adjust adiabatically as time evolves; this implies that the overlap coefficients 
W 2 ** (and also m^*) in (j4~3]) are time-dependent, complicating the extraction of the 
neutron mass shift from the slope of the correlator ratio as a function of time. Even 
with the expanded set of data obtained within the present work, not enough information 
is available to disentangle these time dependences in the most general case. However, 
in one specific situation, which will be argued below to be the relevant one as far as the 
extraction of the neutron electric polarizability is concerned, the slope of R,2{t) indeed 
does yield the neutron mass shift directly; namely, when the Hamiltonian is stationary 
in time. In that case, time dependences in the coefficients W 2 ** (and m") are relegated 
to higher than linear ordei0, and one can indeed equate, up to a minus sign, the slope 
of R.2(t), cf. f l4~3]) . with the mass shift 

Am = m% A A 2 + m* E AE + m EE E 2 . (44) 

Moreover, since a shift in A is equivalent to a shift in time, stationarity of the Hamil- 
tonian in time also implies stationarity in A, i.e., the mass shift (jH|) (and consequently 
the slope of is stationary in A in this particular situation. Thus, in analyzing the 

measured data below, the slope of R2(t), 

will be extracted^ and, for given E, the unique external field will be sought out at which 
5*2 is stationary with respect to A. Since this is then necessarily the point at which the 
Hamiltonian is stationary in time, at that point, then, one can identify 

S 2 = -Am . (46) 



7 As already mentioned at the end of section 12.11 and discussed further in section 14.21 here, an 
additional technical issue arises: While stationarity of the Hamiltonian guarantees that the neutron 
wave function is stationary, there is, in addition, a time dependence contained in the smeared neutron 
sink via eqs. (|28|) .(|29 p . As a result, the overlap between neutron wave function and sink can still contain 
contributions linear in time if one insists on manifest invariance of the neutron sink with respect to 
gauge transformations of the external gauge field, implying the inclusion of (|28p . ([2"9")l in the smeared 
sink construction. On the other hand, if one restricts the calculation to the fixed gauge field (142)) 
and foregoes manifest invariance of the neutron sink with respect to gauge transformations of the 
external field, it is legitimate to use the time-independent smeared sink (|2~T)) alone. In terms of the 
diagrammatic representation of Fig. [TJ this corresponds to discarding all diagrams involving smeared 
sinks other than O- hi the analysis below, both options will be treated, and the final result for the 
neutron electric polarizability will be seen to be uninfluenced by this choice. A way to avoid this issue, 
not explored within the present investigation, would be to use a point neutron sink; such a sink would be 
simultaneously time-independent and invariant under gauge transformations of the external field. On 
the other hand, a point sink would have a small overlap with the true neutron wave function, implying 
a lessened efficiency in the extraction of the neutron ground state signal. 

8 In practice, the average slope over a fixed measurement time interval will be determined in order 
to reduce the statistical uncertainty. 
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To complete the analysis, it is necessary to discuss in more detail the dependence of the 
neutron mass shift Am on the parameters E and A, and, in particular, the relevance of 
stationarity in A. In general, the part of the neutron mass shift which is of second order 
in the external gauge field (j4~2l) can be written as a quadratic form in the parameters E 
and A, cf. flUl) . However, this quadratic form is not yet defined unambiguously and its 
coefficients can consequently not yet all be interpreted as bona fide physical properties of 
the neutron. To see this, consider shifting the entire neutron mass measurement process 
by a time increment i, i.e., the neutron source, which starting with (I35I) has so far been 
assumed to be located at t = 0, shall, for the sake of the following argument, now be 
located at t — i If one concomitantly introduces a shifted time coordinate t' = t — t 
and a shifted 

A = A + Et , (47) 

then, in terms of the shifted quantities, the problem takes a form identical to the original 
one, i.e., one measures the mass shift 

Am = m AA A 2 + m AE AE + m EE E 2 (48) 
= m AA A 2 + (m AE + 2tm AA )AE + (m EE + im AE + 1 2 m AA )E 2 (49) 
= m AA A 2 + m AE AE + rh EE E 2 . (50) 

Thus, in terms of the original definition of E and A, cf. ( 14"2"|) . the shifted measurement 
yields a quadratic form for the mass shift with different coefficients rh AE and m EE 
(whereas the remaining coefficient is invariant, m AA = m AA ). Therefore, the question 
arises how the neutron electric polarizability is to be extracted from the total mass shift 
Am; evidently, polarizability effects enter both the coefficients mf E and rh AE , which 
can be traded off against one another, as demonstrated above. 

As a first step towards disentangling the different effects at play, the immutable 
character of the coefficient m AA should be noted, which allows it to be interpreted as 
an unambiguous property of the neutron. This property moreover is separate from 
the electric polarizability; as verified by explicit calculation below, also at E — one 
obtains the mass shift Am = m AA A 2 , encoding the response of the neutron to distortion 
by the presence of the constant background fielcfi With Am = m A A 2 representing the 
response of a neutron already in the absence of any external electric field, one would 
indeed expect modifications of this response due to a distortion of the neutron by an 
additional electric field to occur only at higher than quadratic order in the external gauge 
field. The representation-independence of m AA , i.e., its independence of the choice off, 
thus seems plausible, and is consistent with the interpretation of the rh AA A 2 term as a 
response separate from the electric polarizability. 

9 Note that the effect of such a constant background field is equivalent to a modification of the 
boundary conditions in the relevant direction, introducing nontrivial Bloch momenta varying with 
quark flavor. 
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On the other hand, in view of (I49p . there is one special representation, i.e., choice 
of i, which seems particularly transparent, namely, the representation in which fh^ 
vanishes, such that! 10 ! 

Am = m* A A 2 + m EE E 2 . (51) 

The form (I5"T1) suggests an interpretation of the data in terms of two, now cleanly 
disentangled effects, namely, the polarizability effect determined by m EE and the effect 
of introducing a constant background field, embodied in the coefficient m^ 4 . Thus, in 
this special representation, the electric polarizability is given by 

a = -2m EE . (52) 

In other words, to isolate the electric polarizability effect from the complete mass shift, 
one simply sets A = in the representation fl5TT) . 

Finally, it is possible to rephrase this prescription for extracting the neutron electric 
polarizability in a manner which is independent of the particular representation, i.e., 
the choice of i. Setting A = in the representation (15TT) is tantamount to evaluating 
the mass shift Am at the extremum in A. However, this way of stating the prescription 
does not rely on that specific representation; after all, in view of (l4~8j) and (!50|) in 
conjunction with (I4"T1) . different representations are related by shifting the value of A, 
and the extremum of Am as a function of A is invariant under such shifts. Thus, one 
can isolate the neutron electric polarizability in any and all representations by seeking 
out the stationary point of the mass shift as a function of A. 

In view of this, and the equivalence of shifts in A with shifts in time, the neutron 
electric polarizability can indeed be extracted by considering the correlator ratio B,2{t), 
cf. (f4"3"j) . specifically for external gauge fields in the vicinity of which the Hamiltonian is 
stationary in time; this validates the arguments presented further above in conjunction 
with eqs. (|H-(|46l). 



3.2 Spatially varying gauge field 

In the case of the gauge field ([3]), 

A = -Ex 3 , (53) 

10 A way to understand how the simplified dependence (|5ip arises is the following: Choosing i such 
as to realize ()51[) shifts the mass shift measurement time interval towards the time t = 0. Now, in view 
of the definition (|42|) . the time t = is special in that the A- and -E-directions in parameter space are, 
in a sense, orthogonal there: At t — 0, a change of E affects only the slope of ^3, but not its value; 
at other times, this is not the case and a change in E also implies an adjustment of the value of A3 
itself, which could be equally effected (or compensated) by a change in A. It is this implicit relation 
between E and A which generates the coupled dependence (|44[) : however, if one measures near t = 0, 
the implicit relation is dissolved and it is natural to obtain the decoupled dependence given by (p)Tj) . 
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one does have a time-independent Hamiltonian, and consequently one can straightfor- 
wardly extract the energy of the neutron ground state from the exponential time decay 
of the neutron two-point function. However, this invariance under translations in time 
comes at the expense of breaking spatial translational invariance. The linear dependence 
of (j55|) on £3 conflicts with the periodic boundary conditions; when traveling through 
the lattice in the 3-direction, as the boundary is traversed, A Q is forced to jump, imply- 
ing a spike in the electric field which is present in addition to the constant electric field 
induced by (1531) . 

Therefore, the neutron is not propagating in a spatially homogeneous background and 
its momentum is not a good quantum numbeJ^l As a result, the ground state energy one 
extracts from the decay of the neutron two-point function contains not only the desired 
mass shift associated with the electric polarizability, but further contributions due to, 
e.g., the effective movement in a spatially varying potential and additional distortions of 
the neutron by the electric field spikes. Within the present investigation, no prescription 
for disentangling the neutron electric polarizability from these other effects with a level of 
cogency comparable to the one discussed in the previous section emerged. Nevertheless, 
the ground state energy obtained below using the external gauge field (j55|) is consistent 
with the polarizability mass shift obtained using the external field (|42|) . suggesting that 
the contamination by the additional effects mentioned above is not dominant. At least 
as far as the neutron's effective propagation in a spatially varying potential is concerned, 
this seems plausible, since the quantum mechanical zero-point energy associated with 
such motion is suppressed by the comparatively large mass of the neutron. 

Note that, although superficially the external fields f|42l) and fl53|) seem quite similar, 
and simply related by an exchange of the temporal with a spatial direction, the physical 
issues arising in the two cases are quite distinct. This is due to the way the mass 
measurement is set up. Up to exponentially suppressed effects, the neutron mass is 
determined by physics within a limited time interval, between neutron source and sink. 
The temporal boundaries, located far in the past or the future of the measurement, 
have a negligible effect on the latter. By contrast, one cannot similarly contain the 
region relevant for the measurement in the spatial directions. In situations with spatial 
translational invariance, by projecting onto a definite momentum, one explicitly weights 
all of space equally during the entire measurement process. Even in the absence of 
spatial translational invariance, it is up to the dynamics to determine whether there 
is a significant probability of finding the neutron near the spatial boundary. Thus, in 
general, the spatial boundary conditions have a crucial influence on the problem. 

In the case of the external field fH2|) discussed in the previous section, this entails 
that shifts of the gauge field A 3 by a constant A have a physical effect, since gauge 

n It is, of course, still legitimate to use a zero-momentum neutron sink, as in (|35[) . since it will 
presumably have a finite overlap with the true neutron ground state wave function. However, that 
wave function itself will not carry a definite momentum. 
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transformations designed to remove such a shift conflict with the spatial boundary con- 
ditions. As a consequence, physics varies locally with time, as discussed extensively 
further above. On the other hand, the neutron mass measurement is insensitive to the 
behavior of the external field at the temporal boundaries. 

In the case of the external field (j55|) . one encounters a largely converse situation: The 
neutron ground state is sensitive to the spatial boundary, at which it encounters spikes 
in the external electric field; on the other hand, as long as one is not in the vicinity 
of the boundary, physics does not vary locally in space. The latter is due to the fact 
that one can indeed remove constant shifts in the field Aq in the time interval relevant 
for the neutron mass measurement using gauge transformations. These transformations 
do of course need to exhibit additional nontrivial structures located far in the past and 
the future of the measurement, but these structures will not influence the measurement. 
This also motivates the fact that no explicit freedom of shifting Aq by a constant is 
included in f l53|) . in contradistinction to fH2|) . Such shifts are not expected to yield new 
physics according to the above argument. 

Comparing the two cases, ultimately (jUJ) can be treated in a more satisfactory 
fashion because the positioning of the neutron source and sink allows one to contain and 
control the breaking of temporal translational invariance introduced by the field (142]) . 
By contrast, in the case of (|53|) . there is no analogous control; the neutron dynamics 
must be allowed to explore space and, in general, the breaking of spatial translational 
invariance will influence the measurement in a nontrivial fashion. 

4 Measurement results 

4.1 Quark wave function renormalization 

To determine the renormalization factor zy in f l2^|) . ff25l) . a measurement of the number 
of valence quarks in the neutron was carried out and subjected to the condition that 
this number equal three. In practice, this is realized by measuring the appropriate 
three-point function, i.e., a diagram of the type 703, with the difference that the lone 
operator insertion is of the form of M 1; cf. ( 1241) , without the weighting by the quark 
electric charge qf, and with a formal external gauge field 



where t is a time between neutron source and sink. Normalizing this by the neutron 
two-point function yields, up to an additional factor i stemming from the Euclidean 
treatment of the time coordinate, the (lattice analogue of the) expectation value of 
/ d 4 x joA in the neutron, where jo denotes the temporal component of the quark current. 
In view of ( 1541) . this reduces to the number of (valence) quarks n = J d 3 x jo present at 
the time t. Fig. [5] displays the plateau obtained measuring n at different insertion times 




(54) 
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Figure 5: (Unrenormalized) number of valence quarks n in the neutron measured for a 
range of insertion times. The neutron source is located at t = and the neutron sink 
at t — 13a. 

t for fixed neutron source and sink. Taking the average of the displayed plateau values, 
one infers 

z v = 1.12 ±0.12 , (55) 

where the uncertainty was obtained using the jackknife method. This measurement of 
zy enters all further measurements below; its uncertainty will be jackknifed into those 
measurements. 

4.2 Neutron mass shift 

Measurements of the neutron mass shift according to the discussion in section [3] were 
taken for the following cases: With the location of the neutron source once again defining 
t = 0, the external gauge field (1421) was studied for the cases to = —10a, to = 
and to — 6a (where a denotes the lattice spacing) at a fixed nonvanishing value of E. 
Furthermore, the case E = at a fixed nonvanishing value of A in (l4"2"j) was investigated. 
The external gauge field (1531) was treated for a fixed nonvanishing value of E, where the 
plane 23 = was taken to define locations maximally distant from the lattice boundary 
in the 3-direction, and was simultaneously used as the location of the neutron source 
(i.e., the smeared quark sources were constructed using an initial position x in (T2T1) 
located in the x 3 = plane). 

Furthermore, as discussed in section 12.31 to suppress fluctuations in disconnected 
diagrams, these external gauge fields were chopped in the time direction, i.e., A% = 
and Aq = for t < —a and t > 14a in the following. Only for —a < t < 14a do A3 and 
A take the forms (1421) and (1531 . respectively. The temporal boundaries of the lattice, 
at which Dirichlet boundary conditions are enforced on the quark fields, are located at 
t = —10a and t = 22a. Stochastic estimation of the disconnected diagrams was based 
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on 120 stochastic sources, as described in section I2.3[ except for the cases t = —10a 
and E = in (T42]) . for which 240 stochastic sources were used. 
Figs. lET-flOl display measurements of the ratio 

= G (p = 0, ( ) ' (56) 

cf. ( I40p . (|43p . for all the aforementioned external gauge fields, in units of the relevant 
external field magnitude. That is, R2 is shown in units of A 2 for the case E — and 
in units of E 2 in the other cases; furthermore, here and in the following, Gaufiian units 
are adopted. Different subsets of diagrams from Fig. [1] contributing to R2 are shown 
in the individual plots (a)-(d) in each case. Figs. l6T-fT0l (a) show only the contributions 
from connected diagrams with lowest-order smearing, i.e., the diagrams 10*. Figs. l6T-fl~0l 
(b) show the result of including all connected diagrams, I**. Note that, in the case of 
the external field (I53p . there are no connected contributions beyond 10*, since smearing 
occurs only in the spatial directions and thus never involves the gauge field component 
Aq. Figs. IBT-TTOl (c) show the result of including all diagrams with lowest-order smearing, 
i.e., the diagrams *0*. Finally, Figs. I6T4T01 (d) show the sum of all diagrams depicted 
in Fig. [TJ Note that, in the SU(3) flavor-symmetric case investigated in this work, the 
only nonvanishing disconnected diagrams are J01 and J03, regardless of the external 
field used. This is due to the fact that the disconnected loop with a linear external field 
insertion is proportional to the sum of the quark charges, and therefore vanishes. Thus, 
Figs. [6HTU1 (c) and Figs. l6T-fT0l (d) contain the same disconnected contributions. 

In comparing Figs. EHHl which display the results obtained using the external gauge 
field (1421) for various t , the different vertical scales should be noted. Compared to the 
case to = —10a, the cases to = and to = 6a exhibit only very small slopes, which 
are determined with relatively small uncertainties, cf . also Table [1] below. The foremost 
observation to be drawn from Figs.[6H8]is that the slope of R2 indeed depends sensitively 
on t , or, equivalently, the constant offset A in the external gauge field fT42|) . as expected. 
This is also corroborated by the E = measurement displayed in Fig. [TU1 As explained 
in section \3A\ these measurements taken together will make it possible to disentangle 
the constant field effect from the electric polarizability. Before proceeding towards this 
central goal, a few further remarks about the data are in order. 

For one, there is a strong cancellation between the diagram 102 and the correspond- 
ing contact term 703, cf. Fig. [TTJ The contact term 703, which originates from expand- 
ing the gauge link variables to second order in the external field, cf. (TTTp . (l2"2"j) . (l2"5"j) . is 
not negligible, as a naive continuum limit might suggest; rather, it contributes to the 
renormalization of the 102 diagram. Presumably, in the continuum limit, it would be 
admissible to disregard diagrams such as 703, at the expense of having to consider a 
strong renormalization of the diagram J02 as its two vertices are permitted to approach 
each other. 
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Figure 6: Contributions by selected subsets of diagrams, as specified in the main text, 
to the ratio i? 2 , as a function of temporal source-sink separation t. All measurements 
are taken at integer times; data are slightly displaced from those times in the figures for 
better readability. These results were obtained using an external field of the form (l42l) . 
i.e., A3 = E(t — to), with to = —10a, where t = corresponds to the neutron source 
location. The electric field E providing the scale is cast in Gaufiian units. 
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Figure 7: Contributions by selected subsets of diagrams, as specified in the main text, 
to the ratio i?2, as a function of temporal source-sink separation t. All measurements 
are taken at integer times; data are slightly displaced from those times in the figures 
for better readability. These results were obtained using an external field of the form 
( l42l) . i.e., A 3 = E(t — t ), with t = 0, where t = corresponds to the neutron source 
location. The electric field E providing the scale is cast in Gaufiian units. 
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Figure 8: Contributions by selected subsets of diagrams, as specified in the main text, 
to the ratio i? 2) as a function of temporal source-sink separation t. All measurements 
are taken at integer times; data are slightly displaced from those times in the figures 
for better readability. These results were obtained using an external field of the form 
( H2l) . i.e., A 3 = E(t — to), with to = 6a, where t — corresponds to the neutron source 
location. The electric field E providing the scale is cast in Gaufiian units. 



0.5 



0.25 



CO 

CM 

DC 



-0.25 



-0.5 '— 1 



"i 1 1 r 

subset (a) i — •— 



• » a 



i 1 r 



"■'if 



0.25 



CO 

C\J 

LU 



DC 







6 

t/a 



10 



12 




-0.25 



-0.5 L - 1 



Figure 9: Contributions by selected subsets of diagrams, as specified in the main text, 
to the ratio R2, as a function of temporal source-sink separation t. These results were 
obtained using an external field of the form (|53l) . i.e., A Q = —Ex^, for this background, 
there are no smearing contributions beyond zeroth order in the external field, i.e., case 
(b) is identical to case (a), and case (d) is identical to case (c). The electric field E 
providing the scale is cast in Gaufiian units. 
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Figure 10: Contributions by selected subsets of diagrams, as specified in the main text, 
to the ratio i?2, as a function of temporal source-sink separation t. All measurements 
are taken at integer times; data are slightly displaced from those times in the figures for 
better readability. These results were obtained using an external field of the form fj42l) 
with E = 0, i.e., A 3 = A. The constant background gauge field A providing the scale is 
cast in Gaufiian units. 
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Figure 11: Contributions by diagrams 102 and 703, as well as their sum, to the ratio 
i?2, as a function of temporal source-sink separation t. All measurements are taken 
at integer times; data are slightly displaced from those times in the figure for better 
readability. These results were obtained using an external field of the form (H2"l) . i.e., 
A 3 = E(t — to), with t = —10a, where t = corresponds to the neutron source location. 
The electric field E providing the scale is cast in Gaufiian units. 
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Secondly, one can furthermore observe from Figs. IqT-TTUI that the disconnected dia- 
grams consistently tend to give a negative contribution to the slope of the ratio 
for an external field of the form (T4*2l . cf. also Table [1] below. For the external field fl53|) . 
the contribution is very slightly positive. However, it should be emphasized that the 
contributions are in no case large enough to be significant compared to the statistical 
uncertainty. 

Thirdly, as noted in section [3TT| the higher-order (in the external field) sink smearing 
diagrams contained in Figs. [6HH] (b),(d) in general contribute additional linear time 
dependences to the ratio R2(t) beyond the ones associated with the mass shift of the 
neutron in the external field. Indeed, a substantial difference in slope can be seerJ^I 
comparing Fig. (a) and (b), or also Fig. (c) and (d). To isolate the slope due to 
the mass shift itself, one should evaluate Figs. [6HH] (c); this comes at the expense of 
foregoing a form of the neutron sink (and source) which is manifestly invariant under 
gauge transformations of the external field. Of course, it is not imperative to use a 
manifestly invariant form; all that is implied by not doing so is that the neutron mass 
shift, a gauge- invariant quantity, has been evaluated in a specific gauge for the external 
electromagnetic field. Note, moreover, that this point is largely moot at any rate, since 
the final result for the neutron electric polarizability will be seen to not be affected 
significantly by the differences in intermediate data introduced by the aforementioned 
sink smearing effects. 

Returning to the main objective, extracting the neutron mass shift from the slope of 
the ratio R2{t), Table [JJ lists the slopes extracted from the data displayed in Figs. IdT-TTUI 
These slopes were obtained by performing least-square fits of linear functions in t to 
the R2{t) data for a range of t; the uncertainties were obtained by jackknife analysis. 
The time range used was 4a < t < 10a; this choice minimizes statistical uncertainty (by 
using as large a time range as possible) while still allowing for a good least-square fit 
by a linear form as well as a jackknife bias estimate small compared to the statistical 
uncertainty^!. 

According to the discussion in section 13. 1[ to extract the electric polarizability from 
the data in the first three lines of Table [TJ these data should be viewed as defining a 
parabola in t Q , and the extremum in i should be sought out (note that this is equivalent 
to viewing the data as defining a parabola in A, since A = —Eto and E is constant). 
At the extremum, one can then identify the slope S2 with the (negative) neutron mass 

12 By considering the diagrams individually, one can indeed verify that, as expected, the difference 
arises specifically due to sink smearing contributions, and not source smearing contributions. 

13 In the case of the external field (|42| with to = and to = 6a, a more restricted time range improves 
the linear fit, as is apparent, e.g., from Fig. [8l however, as already noted above, these two instances are 
determined with very little uncertainty (regardless of the time range used) compared to to = —10a. As 
a result, it is the latter case which dictates the choice of time range, which then was adopted for all 
cases for consistency. 
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(a) 


(b) 


(c) 


(d) 


m,t = -10a 


0.46(18) 


0.83(24) 


0.26(26) 


0.63(30) 


S 2 /{a»E*) 
(42),t = 


0.000(16) 


0.096(38) 


-0.033(43) 


0.063(54) 


S 2 /{a*E*) 
(H2D, t = 6a 


-0.017(3) 


-0.047(7) 


-0.037(25) 


-0.067(27) 


-Am/{a 6 E A ) 

(EHD 


-0.027(4) 


-0.027(4) 


-0.019(41) 


-0.019(41) 


-Am/(aA 2 ) 
dH, £ = 


0.0025(8) 


0.0028(8) 


0.0017(9) 


0.0020(10) 



Table 1: Slope S 2 , cf. (|45p . for different external electromagnetic fields, in the appropriate 
external field units and units of the lattice spacing a. Cases (a)-(d) correspond to the 
different subsets of diagrams included in the corresponding Figs. [BlfTUl (a)-(d), cf. main 
text. In the case of the background field (1531) . as well as for ( 1421) with E = 0, the slope 
S 2 can be directly identified with the negative mass shift, —Am; hence the labeling 
of the last two lines. By contrast, the data in the first three lines must be processed 
further to locate the stationary point as a function of to, at which then S 2 = —Am can 
be identified, cf. main text and Table [2j 



shift, -Am. Fitting the form (cf. pE]) - fl4"B]) ) 

S 2 /E 2 = -(m AA A 2 + m AE AE + m EE E 2 )/E 2 

= - (m AA t 2 - m A % + m EE ) (57) 

to the data in Tabled] yields parabolas with the extrema and curvatures listed in Table [2j 
The uncertainties quoted in Table [2] were again obtained using the jackknife method. 
The quadratic coefficient m AA extracted in this way agrees well with the E = valued 
listed in the last line of Table [H providing an independent measurement corroborating 
the interpretation of the data advanced in section 13.11 

The central result of this work, however, is the value of the electric polarizability of 
the neutron, the full value of which is obtained by multiplying the mass shift quoted 
in column (c) of Table [2] by a factor of 2, cf. (|4"4"|) . (1521) . In physical units, obtained by 
inserting a = 0.124fm, one has 

a = —2Am/E 2 = (-2.0 ± 0.9) ■ 10~ 4 fm 3 . (58) 
14 The last lines of Tables [T] and [5] can be directly compared, since for E = 0, one has Am/ A 2 = m^ A , 

cf. 
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(a) 


(b) 


(c) 


(d) 


-Am/(a 3 E 2 ) 
(extremum) 


-0.034(6) 


-0.049(10) 


-0.052(24) 


-0.076(63) 


—vn^ /a 


0.0027(9) 


0.0031(9) 


0.0018(16) 


0.0022(16) 



Table 2: (Negative) mass shifts at the extrema of the parabolas defined by the data 
in Table [T] through the form (|57p . as well as the coefficient characterizing the 

curvatures of the parabolas. 

Note that the additional smearing contributions entering the result in column (d) indeed 
do not significantly alter this result. The result (1581) is corroborated by the measurement 
using the external field (153"|) , quoted in the fourth line of Table [TJ translated into physical 
units, that measurement would imply a polarizability of 

a = (-0.7 ± 1.6) • 10" 4 fm 3 . (59) 

As discussed in section I3.2[ the result fl59|) contains systematic uncertainties (not in- 
cluded in the quoted statistical error) stemming from the fact that the mass shift mea- 
sured in this case is contaminated by the quantum mechanical zero-point motion of the 
neutron and distortions of its internal wave function due to superfluous spikes in the 
external electric field. Thus, the result (1581 is expected to be more trustworthy than the 
result (159"]) . Nevertheless, the difference between the two measurements does not turn 
out to be significant; the aforementioned contaminations do not appear to represent 
appreciable effects. 

Compared to the experimental value reported by the Particle Data Group [25] 

a = (11.6 ± 1.5) • 10~ 4 fm 3 (60) 

the result (|58|) suggests a strong variation of the electric polarizability of the neutron 
with the pion mass. Indeed, Chiral Effective Theory calls for such a variation [10,26-29], 
dominated by a l/m^ dependence at low pion masses. In the "Small Scale Expansion" 
approach [27,28], which systematically extends leading-one- loop Heavy Baryon Chiral 
Perturbation Theory by including explicit A degrees of freedom, the electric polariz- 
ability of the neutron decreases by an order of magnitude as one varies the pion mass 
from the physical point up to around 400 MeV. Qualitatively, a change of sign of the 
polarizability at even higher pion masses, as implied by f[55|) . does not seem implau- 
sibld^l, although it should be stressed that a pion mass of 759 MeV, corresponding to 
the dynamical quark ensemble used in the present work, is certainly far beyond the 

15 Of course, for very large masses, i.e., in the nonrelativistic limit, the polarizability cannot be 
negative due to the general properties of second-order perturbation theory. 
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regime in which Chiral Effective Theory can be applied reliably. Lattice calculations at 
lower pion masses are needed in order to achieve a quantitative connection with Chiral 
Effective Theory. 

On the other hand, the result obtained in the present work at first sight appears 
to be at odds with previous lattice measurements [5,7]. Those studies yield a neutron 
electricpolarizability which is consistent with f[6"Ul) over a wide range of (valence) quark 
masses^j, including the quark mass used in the present work. The question arises how 
such a weak variation with the quark mass can be reconciled with the result obtained in 
the present work (and also with the expectation coming from Chiral Effective Theory). 

Apart from the use of the quenched approximation, the main differences between [5,7] 
and the present treatment are that, on the one hand, [5, 7] work with an external field 
corresponding to a particular value of to in (j42l) ; on the other hand, this external field 
is introduced into the lattice link variables in linearized form, i.e., the vertex insertion 
M 2 , cf. (125]) . is not included. As discussed further above, the insertion M 2 provides con- 
tact terms which renormalize propagators with two Mi insertions, leading to substantial 
cancellations; however, most importantly, the measured mass shift depends sensitively 
on the parameter to characterizing the external field. It is instructive to reevaluate the 
data gathered in the present work such that diagrams generated by M 2 insertions are 
excluded, and at a value of to corresponding to the one used in [5,7]. Two measurements 
are provided by [5,7]. One uses Wilson fermions with a distance of 1.7 fm between t and 
the mass shift measurement, yielding a polarizabilit}0 of a = (9.8 ± 1.2) • KT 4 fm 3 for 
pion masses comparable to the one at which the present work was performed; the other 
uses clover fermions with a distance of 1.53 fm between to and the mass shift measure- 
ment, yielding a polarizability of a = (13.9 ± 0.8) • 10~ 4 fm 3 at comparable pion masses. 
Taking into account that the mass shift measurement in the present work is centered 
around t = 7a = 0.87 fm, the values of t corresponding to the two aforementioned cases 
are t = —6.7a and t = —5.3a, respectively. Constructing the parabola floTl) defined 
by excluding diagrams generated by M 2 insertions from the set of diagrams comprising 
case (c) above, and evaluating it at those values of to yields 

a(t = -6.7a) = (20 ± 11) • 10" 4 fm 3 (61) 
a(t = -5.3a) = (15 ± 8) • 10~ 4 fm 3 (62) 

which, particularly in the latter case, is in quite good agreement^ with the results 
of [5,7]. Thus, at the level of the raw numerical measurement, the present work in 
fact corroborates the results obtained in [5,7]; at the same time, it is now clear that 
such a measurement at a single fixed to in general contains two separate effects, i.e., 

16 The studies reported in [5,7] employ the quenched approximation. 
17 These estimates were obtained by linearly interpolating results quoted in [7] in m^. 
18 It should be noted that these comparisons depend sensitively on the determinations of the lattice 
spacings in the different calculations. 
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the electric polarizability itself and the effect of subjecting the neutron to a constant 
external field. To disentangle the former from the latter using such measurements, it is 
necessary to combine the data obtained using a variety of t . This suggests that it would 
be worthwhile to supplement the measurements already performed in [5, 7] by further 
analogous measurements at other values of to in order to obtain a more comprehensive 
picture of the physical effects engendered by the introduction of the external field and 
thus be able to isolate the different effects from one another. 

5 Summary and outlook 

The investigation reported here represents a first exploration of the neutron electric 
polarizability in the context of lattice QCD with dynamical quarks. Its main thrust lay 
in clarifying conceptual questions within the framework of the background field method 
and assessing the feasibility of numerical computations on that basis, using a SU(3) 
flavor- symmetric ensemble as a test case. 

Two central issues needed to be addressed to arrive at a cogent calculational scheme. 
On the one hand, the presence of dynamical quarks dictates the use of four-point func- 
tion methods, introducing, in particular, the need to evaluate disconnected diagrams. 
These contributions, which were included in the numerical calculations carried out in 
this work via stochastic estimation, significantly increase the computational expense of 
the measurement. Nevertheless, the feasibility of carrying out such measurements was 
demonstrated for the SU (3) flavor-symmetric case; the cost of progress towards lighter 
quark masses does not seem prohibitive, but such an endeavor will require a, by current 
standards, significant commitment of computational resources. 

On the other hand, a strong emphasis was placed within the present work on the 
physical consequences of shifting the external electromagnetic field by a constant on 
a finite lattice. While such shifts merely correspond to gauge transformations in infi- 
nite space, on a finite lattice, they influence the physical spectrum and thus mask the 
mass shift due to the electric polarizability itself. On lattices of a practical size, this 
effect has a dominating influence on the neutron mass shift, from which one aims to 
extract the electrical polarizability. To disentangle the two effects, measurements using 
a variety of external fields which are shifted with respect to one another are necessary 
(further impacting computational cost). It should be noted that this issue also affects 
investigations carried out in the quenched approximation, such as reported in [5,7]. The 
present investigation, complementing that effort, suggests additional measurements to 
supplement the ones already carried out, in order to gain a comprehensive picture of the 
effects playing a role. It is hoped that the results obtained here will provide motivation 
and useful input for an expanded measurement program in this direction. 

Looking forward, besides the obvious need to progress towards lighter quark masses, 
it would be interesting to study other hadrons, especially with a view towards measuring 
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polarizability combinations in which disconnected diagrams at least partially cancel. 
Such combinations could be calculated with higher accuracy at lower cost. However, 
a potential obstacle to this which should be kept in mind is the following: Typically 
hadrons of differing electric charge would be involved, and, a priori, it is not clear that 
measurements using the same external electromagnetic field are appropriate in each 
case for the purpose of isolating the electric polarizability. On the other hand, results 
obtained in different external fields cannot be combined straightforwardly to cancel 
disconnected contributions. 

One possibility of avoiding such difficulties lies in using alternate methods of access- 
ing polarizabilities, e.g., via density- density correlation functions [30]. That approach 
would circumvent the necessity of explicitly introducing an external electromagnetic 
field. Density-density correlation functions at unequal times can be used to extract 
hadron polarizabilities, specifically by measuring the second moment (with respect to 
spatial separation) of the correlation function for a range of relative times and inte- 
grating over the latter. When calculating polarizabilities of hadrons in this manner, at 
least partial cancellations of disconnected diagrams can be achieved straightforwardly 
by forming the proper isovector combinations. Care must be taken to restrict the hadron 
momentum to the nonrelativistic regime, in order to exclude relativistic effects which 
complicate the interpretation of the density-density correlation function and the extrac- 
tion of the polarizability. Also, density- density correlation functions generally fall off less 
rapidly than standard hadron wave functions (a doubling of the extent being typical); 
this has motivated the development of periodic image correction methods [30] which are 
expected to prove helpful in this context. 
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and also in [2,5,7] should be analytically continued to imaginary values. This further 
step, which was not carried out in [2, 5, 7] and in the present treatment, implies that 
the results (158]) and fl59l) receive an additional overall minus sign. In view of the small 
magnitude of fl58l) and fl59l) . this does not decisively impact the further conclusions 
drawn in the present work. However, it will need to be taken into account, and revisited 
in more detail, in work going forward aiming at progress towards lighter pion masses. 
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